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ABSTRACT 

We report the detection of off-pulse emission from two long period pulsars PSR B0525+21 
& PSR B2045-16 using the Giant Metrewave Radio Telescope (GMRT). We recorded high time 
resolution interferometric data at 325 MHz and gated the data offline to separate the on- and 
off-pulse sections of the pulsar. On imaging the two temporal sections separately, we detected a 
point source in the off-pulse images of both the pulsars coincident with the pulsar locations in 
the on-pulse images. The long periods (3.75 and I .96 s respectively) and low spin-down energies 
(E ~ fO 31 erg s _1 ) of the two pulsars argue against a pulsar wind nebula (PWN) origin for the 
off-pulse emission. The derived estimates of the ISM particle density required to drive a PWN 
do not support such an interpretation. A magnetospheric origin for the off-pulse emission raises 
questions regarding the location of the emission region, which is an important input into pulsar 
emission models. 

Subject headings: pulsars: general — pulsars: individual (B0525+21, B2045-16) 



1. Introduction 

We report the detection of off-pulse emission 
from two long period pulsars PSR B0525+2I and 
PSR B2045-I6 using the interferometric mode of 
the Giant Metrewave Radio Telescope (GMRT) in 
the 325 MHz frequency band. The "main-pulse" 
(or the equivalent term "on-pulse" used in this pa- 
per) is the emission from within the polar-cap re- 



1 Part of the work was done while at NAICenter, Arecibo 
Observatory 



gion of a pulsar; and the off-pulse is the emission 
outside this main-pulse. Off-pulse emission from 
pulsars has been a subject of interest since the 
discovery of pulsars four decades ago. 

According to the rotating vector model (RVM) 
proposed by Radhakrishnan & Cooke (1969) the 
observed pulses are due to relativistically beamed 
radiation along the open dipolar field lines. The 
plane of linear polarization traces the magnetic 
field line associated with the emission at every 
instant. Several statistical studies have revealed 
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that pulsar radio emission is in the form of a cir- 
cular beam arising from the polar cap (i.e. the 
main-pulse) with an opening angle p, and from a 
height above the magnetic pole equal to 1-2 per 
cent of the radius of the light cylinder. The width 
of the observed pulse is a geometrical function of p 
(which itself is a function of the emission height), 
the pulsar period, the angle between the rotation 
axis and the magnetic axis (a) and the angle be- 
tween the magnetic axis and the line of sight to 
the observer (0). The angles a and j3 can be es- 
timated by fitting the RVM to the swing of the 
polarisation position angle (PPA) across the pulse 
([Everett fc Weisberg 200ll iMitra fc Li 20041 ). 

The observed main-pulse covers 5-10 per cent 
of the pulse period for 90 per cent of the pul- 
sars. In several pulsars this main-pulse consists 
of multiple components including low-level bridge 
emission between such components (Rathnasree & 
Rankin 1995). The pulse width scales as l/sin(a), 
from the geometry of the emission model; in the 
rare case of a pulsar with closely aligned mag- 
netic and rotation axes (i.e. a ~ 0° or 180°) the 
pulse width can be as high as 100 per cent. Sev- 
eral pulsars have low level but detectable bridge 
emission between two widely spaced pulse com- 
ponents, making them candidates for an aligned 
rotator geometry. In rare cases when a is close 
to 90°, emission from the opposite pole (i.e. the 
inter-pulse) may be observed at 180°phase from 
the main-pulse. In all these instances the observed 
emission is believed to lie within the polar cap (i.e. 
main-pulse) as described earlier. 

In a few pulsars, low level emission components 
known as pre-/post-cursors (PPC) have been ob- 
served outside the main-pulse (e.g. see Mitra & 
Rankin 2010 for a discussion). In these pulsars, 
all geometrical evidence indicate that the main 
pulse is consistent with emission from open field 
lines. The PPC components appear highly polar- 
ized and are far from the main-pulse. The discov- 
ery of the PPC component about 60°away from 
the main pulse in PSR B0943+10 is particularly 
interesting, as the line of sig ht almost grazes th e 
emission cone for this pulsar ( Backus et al. 201oh . 
Hence the PPC emission originates either from 
a much larger height, where due to spreading of 
dipolar field lines the PPC component can lie far 
away from the main-pulse, or from the regions of 
closed field lines. Some pulsars occasionally emit 



giant pulses which are believed to arise close to 
the light cylinder rather than the polar cap. In- 
terestingly PPC pulsars have low spin-down lu- 
minosities (E < 10 34 erg s _1 ), whereas the giant 
pulses arise from pulsars with E > 10 34 erg s _1 . 
The PPC components and giant pulses are poten- 
tial sources of magnetospheric off-pulse emission in 
pulsars. They challenge the conventional wisdom 
of pulsar radio emission arising only from open 
magnetic field lines close to the pulsar polar cap 
and raise questions about the origin of these emis- 
sion components. 

The other possible source of off-pulse emission 
is the interstellar medium (ISM) around the pul- 
sar. The pulsar loses most of its rotational en- 
ergy in the form of a relativistic wind which, when 
confined by the surroundings, may form a pulsar 
wind nebula (PWN). Several types of PWNe are 
formed depending on the confinement mechanism. 
Young and energetic pulsars are often located in 
their associated supernova remnant (SNR). The 
pulsar wind streaming into the ambient medium 
produces standing shocks resulting in a plerionic 
PWN (like the one observed in the Crab nebula). 
In older pulsars, where the surrounding SNR is 
likely to be dissipated, the relativistic particles 
may interact with the ISM magnetic field and radi- 
ate, creating ghost remnants (proposed by Bland- 
ford et al. 1973, but not observed till date). Fi- 
nally, a pulsar moving through the ISM with su- 
personic speed can produce a bow-shock nebula, 
where ram pressure balance is established between 
the pulsar wind and the ambient medium. Several 
such bow shock nebulae have been detected in Ha, 
and some in the radio as plerionic bow shocks, but 
never in both (see Chatterjee & Cordes 2002 for a 
detailed study). 

Several searches for PWNe have been con- 
ducted in the past with varying sensitivities and 
resolutions. However only in about 10 young pul- 
sars have radio emission been detected outside 
the main-pulse and all of these are believed to 
be associated with PWNe. In all such cases the 
associated pulsars are young (10 3 -10 5 yr), have 
high spin-down luminosities (E > 10 35 erg s _1 ) 

and all except one are associated with SNRs 

(|Gaensler et al. 1998tlGaensler et al. 2000tlstappers et al. 19991) . 
It is believed that the pulsar wind and the envi- 
ronment change as pulsars slow down and age, 
making them less likely to harbour PWNe. 
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The present study was designed to detect radio 
emission away from the polar caps and not associ- 
ated with PWNe. Therefore we targeted old and 
less energetic pulsars unlikely to harbour PWNe. 
Their profiles did not show any components out- 
side the main-pulse. The targets were neither 
aligned rotators nor were inter-pulsars. We at- 
tempted to detect emission within a pulse phase 
range of approximately 80-250° from the peak in 
the main-pulse. 

We describe the gated interferometric observa- 
tions and the data analysis in section 2, the results 
in section 3, the additional tests we carried out to 
confirm that the off-pulse emission were not in- 
strumental artifacts in section 4, and discuss the 
implications of the detections in section 5. 

2. Observations & Analysis 

2.1. Target Pulsars 

We chose two pulsars, B0525+21 and B2045- 
16, whose pulse profiles did not show any fea tures 
outside the main-pulse (jGould fc Lvne 19 98) and 
were stronger than 50 mJy to be able to put a 
non-detection upper limit of at most 5 per cent of 
the pulsed flux. The properties of the targets are 
listed in table [TJ We targeted long period pulsars 
because: 

1 . PWNe are so far only known from short-period, 
energetic pulsars in the vicinity of SNR. A long pe- 
riod selects against both high energy pulsars and 
SNR association, and hence PWNe. 

2. We wanted at least 8 bins (of 131 or 262 ms 
each) across the pulse period to cleanly separate 
the off-pulse and the on-pulse regions. 

The temporal broadening across the 16 MHz 
bandwidth due to dispersion by the ISM was 183.3 
ms for PSR B0525+21 (binwidth 262 ms, gate 
width 1.31 s), and 47.8 ms for PSR B2045-16 (bin- 
width 131 ms, gate width 0.655 s). Therefore, we 
did not have to dedisperse the signal. 

2.2. Interferometric imaging of the pulsars 
with the GMRT 

Interferometric observations are better than 
single dish measurements for this study for sev- 
eral reasons: 

1 . An imaging interferometer is only insensitive to 
constant flux density background along the spa- 



tial axes while a standard pulsar receiver is usu- 
ally insensitive to the constant background along 
the time axis; the detection of off-pulse emission 
is essentially an attempt to find such a constant 
background along the time axis. 

2. Self-calibration of interferometric data can cor- 
rect instrumental and atmospheric gain fluctua- 
tions on very short time scales. The corrections 
arc determined by the flux densities of the con- 
stant and bright background sources in the field 
and hence would not be affected by the pulse vari- 
ation of the relatively weak pulsar. 

3. The higher spatial resolution of an interfer- 
ometer greatly reduces the coincidence of unre- 
lated sources; thereby reducing the probability 
that the off-pulse emission is from an unrelated 
source within the synthesized beam. 

We imaged the pulsars with the Giant Me- 
trewave Radio Telescope (GMRT), an aperture- 
synth esis radio interferom eter located near Pune, 
India (jSwarup et al. 19911 ). The 30 antennas of 45 
m diameter provide a maximum baseline of 27 km 
and can be operated at 6 frequency bands between 
50 and 1450 MHz. We observed the pulsars at 
325 MHz with a 16 MHz bandwidth split into 128 
channels. The frequency was chosen for its opti- 
mal combination of resolution (10 ") and sensitiv- 
ity (few 100 fiJy in 4 hours of observing) . The high 
frequency-resolution was useful in flagging narrow- 
band radio-frequency interference (RFI) and to 
avoid bandwidth smearing. The shortest integra- 
tion output by the hardware correlator is 0.131 s, 
though this data is usually averaged and median 
filtered to provide a standard data output integra- 
tion of 16.77s. 

The two pulsars were first observed on 19 May, 
2009 with 0.262s integration to keep the data out- 
put rate and size manageable. This time resolu- 
tion (of 7 bins) was found to be inadequate in the 
case of PSR B2045-16 while the data was of poor 
quality (> 10 missing antennas) for both targets. 
We reobserved the sources on 19 January, 2010, 
with integrations of 0.131s for PSR B2045-16 (pe- 
riod 1.9616s), and 0.262s for PSR B0525+21 (pe- 
riod 3.7455s). This allowed us to divide both the 
pulsar periods into 14 time bins. 

2.3. Data Analysis 

Several earlier studies used the technique of 
online gating to separate the on- and off-pulse 
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Table 1: Properties of Pulsars selected for off-pulse studies. 



Pulsar 




period 


DM 


dist 




E 


Vtrans 


Flux 32 5 






s 


cm~ 3 pc 


kpc 


year 


erg s _1 


km 


mjy 


B0525- 


1-21 


3.7455 


50.937 


2.28 


1.48x10" 


3.0xlO ai 


229 


80.5 


B2045- 


-16 


1.9616 


11.456 


0.95 


2.84xl0 6 


5.7xl0 31 


511 


169.3 



Period, dispersion measure and proper motion (for Vtrans) are from (|Hobbs et al. 200^); the oth ers are from 
(|Tavlor et al. 1993h . FIUX325 is the expected flux at 325 MHz calculated from values in l Lorimer 1995h . 



regions using the timing information of pul- 
sars, and acc umulated time-averaged on- and off- 
pulse data (IStrom fc Van Someren Greve 1 990; 



Stappers et al. 19991 ) . This resulted in relatively 
smaller size of data files and hence shorter com- 
putation times which was a major consideration 
1-2 decades ago. The disadvantage was the impos- 
sibility of rectifying any errors during the gating 
process (see results and discussions in Gaensler et 
al. 2000). We recorded the data at sufficiently 
high resolution and applied offline gating. 

2.3.1. Folding & Gating 

In the absence of accurate absolute time refer- 
ence in the interferometric mode we used offline 
gating to isolate the off-pulse region of the pulsar 
period. The 30 antennas of GMRT provide 30 self, 
and 435 cross visibilities at each integration inter- 
val (0.131s or 0.262s). We folded the self-data to 
obtain the pulsar profile and identify the on-and 
off-pulse regions. This was used to gate the cross- 
antenna visibilities into on- and off-pulse data sets. 

The temporal and frequency gain variations in 
the self-data were corrected by normalising the 
instantaneous values by the local median. The 
timescale for the temporal median was 30-50 times 
the pulsar period to ensure that the pulse varia- 
tion was retained. The data from all frequency 
channels and antennas were averaged with robust 
sigma-clipping and processed using the standard 
folding algorithm for pulsars (Hankins & Rickett 
1975) to determine the pulse profile. The pulsar 
profiles and the off- and on-pulse gated are shown 
in Figure [TJ 

The off- and on-pulse gates of PSR B0525+21 
were centred on phases 0.24 and 0.64, respectively. 
The 5 bins nearest to each were averaged to con- 
struct the off- and on-pulse data sets. The corre- 



sponding phases for PSR B2045-16 were 0.21 and 
0.67, respectively. The folding and gating proce- 
dures were carried out using software developed 
for this purpose. 

2.3.2. Imaging 

The on- and off-pulse data were cleaned of radio 
freque ncy interferenc e (RFI) using the RfiX algo- 
rithm ( Athreva 20091 ). They were separately cali- 
brated, flagged and imaged in a standard manner 
using the NRAO AIPS package. The flux-density 
scale of the images was determined from observa- 
tions of the calibrator 3C48 and the measurements 
of Baars et al.(1977) with the latest corrections 
of Perley et al (1999; in AIPS). The data sets 
were taken through several rounds of phase self- 
calibration followed by a final round of amplitude 
and phase self-calibration. 

PSR B0525+21: Observations of the calibra- 
tor 0521+166 were interspersed with the target 
to correct for the amplitude and phase gain fluc- 
tuations. The presence of the extremely strong 
and extended crab nebula (flux > 1000 Jy and 
angular size ~ 10') at the edge of the primary 
beam (1.5° away) resulted in enhanced noise and 
strong ripples in the initial image. Therefore, sub- 
sequent analyses were carried out with a lower UV 
cut-off at 1.5 kA which reduced the artifacts and 
noise. This eliminated all structures larger than 
2.3' from the image; this did not affect our pur- 
pose as our target source was expected to be much 
smaller than 2'. 

PSR B2045 16: The calibrator 2137-207 was 
observed for amplitude and phase calibration. A 
pointing error resulted in the pulsar being located 
35' away from the field centre (60 per cent gain 
level of the primary beam). The primary beam 
correction was applied to the image to obtain the 
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B0525+21 408 MHz (Gould & Lyne, 1998) 
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1. — Pulse profiles for our targets from standard pulsar mode high time resolution observations (top 
and our folded interferometric self data (bottom row). 
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correct flux density for the pulsar. 



3. Results 



We detected off-pulse emission in both the pul- 
sars, coincident with the location in the on-pulse 
image. The results are summarized in tabic [2] 
The on- and the off-pulse emissions appeared to 
be unresolved for both the pulsars. 

We compared the flux densities of sources in the 
field of view to check the flux-density calibration. 
The ratios of on- and off-image flux densities of 
sources stronger than 10 mJy are plotted in figure 
2J There were 43 and 93 comparison sources in 
the fields of B0525+21 and B2045-16, respectively. 
The scatter of values around unity, especially for 
strong sources (which have smaller fractional flux 
errors) confirmed the similarity of the flux scale of 
the on- and off-images for both the pulsars. 

PSR B0525+21: The pulsar was detected in 
the off-pulse image at the 8.6a level (flux-density 
3.9±0.5 mJy, rms noise 0.45 mJy/beam). The flux 
density of the pulsar in the on-pulse image was 
80.2±5.8 mJy (rms noise 0.55 mJy/beam). The 
resolution element at the distance of the pulsar 
corresponds to 0.09 pc. The data from 19 May, 
2009, yielded an off-pulse detection of 3.6 mJy. 

PSR B2045 16: The pulsar was detected in 
the off-pulse image at the 6.6a level (flux density 
4.3±1.1 mJy; rms noise 0.65 mJy/beam). The flux 
density of the pulsar in the on-pulse image was 
305.2±21.9 mJy (rms noise 2.1 mJy/beam). The 
resolution element at the distance of the pulsar 
corresponds to 0.04 pc. 

The flux-densities averaged over the entire 
pulsar period is 30.0±2.1 for B0525+21 and 
110.5±7.9 for B2045-16. These values are smaller 
than expected (Table [1]) by factors of 1.5 and 2.7, 
which can be caused by effects such as refractive 
interstellar scintillations (jStinebring fc Condon 199ph . 
We confirmed the correctness of our flux scale by 
comparing the flux densities of the phase calibra- 
tors to their known values (less than 1 per cent) 
and by comparing the flux densities of 3-5 strong 
sources in each field to the interpolated values 
from NVSS (1.4 GHz) and VLSS (74 MHz). 

The quoted errors on the flux densities were 
obtained by adding in quadrature the image rms, 
calibration errors and the error on the value of the 



PSR B0525+21 On pulse 



22 01 15 

00 

00 45 

f 30 

i 

i 

'. 15 



w 21 59 45 
30 
15 

00 - 
58 45 



i =j i i r 




05 28 58 56 54 52 50 48 

RIGHT ASCENSION (J2000) 
Levs = 0.55 mJy * (-5, -3, 3, 4, 5, 10, 20, 40, 60, 80, 100, 120) 

flux = 80.18 mJy, rms = 0.54 mJy/beam 



PSRB0525+21 Off pulse 



22 01 15 
00 
00 45 



< 00 

o 

g 21 59 45 
30 



58 45 



T 



56 



46 



05 28 58 56 54 52 50 

RIGHT ASCENSION (J2000) 
Levs = 0.45 mJy * (-5, -3, 3, 4, 5, 6, 7) 
flux = 3.87 mJy, rms = 0.45 mJy/beam 



Fig. 2. — Contour plots of the on- and off-pulse 
images of B0525+21 showing the pulsar detection. 
The synthesized beam is shown in box. 



primary flux-density calibrator (jBaars et al. 19771 ) 
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Table 2: Summary of on- and off-pulse observations and results. 



Pulsar 




frequency 


time 


beam 


size 


On rms 


On Flux 


Off rms 


Off Flux 


Int. Flux 






(MHz) ' 


(min) 


(") 


(PC) 


(mJyb- 1 ) 


(mJyfe- 1 ) 


(mJyb- 1 ) 


(mJyfc- 1 ) 


(mJyfe- 1 ) 


B0525- 


1-21 


332.9 


160 


9.5x6.5 


0.088 


0.55 


80.2±5.8 


0.45 


3.9±0.5 


30.0±2.1 


B2045- 


-16 


317.1 


180 


11.9x7.2 


0.042 


2.10 


305±22 


0.65 


4.3±1.1 


110.5±7.9 
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Fig. 4. — The ratio of the measured flux-densities 
of point sources (flux density > 10 mjy) in the on- 
and off-pulse images. The scatter around unity 
confirmed that the flux-scale was the same for 
both. (B0525+21: top; B2045-16: bottom) 



Fig. 3. — Contour plots of the on- and off-pulse 
images of B2045-16 showing the pulsar detection. 
The synthesized beam is shown in box. 
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The signal-to-n oise ratio in the publ ished pro- 
files of the pulsars ( Gould fc Lvne 19981 ) are insuf- 
ficient to conclude whether the off-pulse emission 
detected here is a pedestal (throughout the period) 
or a pulse confined to the off-pulse gate. 

4. Genuineness of the Off-Pulse De- 
tection 

Firstly, we detected the off-pulse emission in 
two data sets of B0525+21 observed 8 months 
apart, and within the noise the two measurements 
were identical. We can safely rule out the detec- 
tions being chance occurance of noise peaks. Sec- 
ondly, B2045-16 was positioned almost at the half- 
power point of the primary beam while B0525+21 
was close to the field centre. The detection in both 
cases suggests that they are not the spurious struc- 
tures occasionally seen at the centre of the field. 
We discuss below several effects which could result 
in spurious sources in the off-pulse data at the lo- 
cation of the pulsar; and demonstrate that they 
are unlikely to be responsible for the same. 

4.1. Positional coincidence of an unrelated 



4.2. Error in Time stamp 

This scenario requires that approximately 1-4 
per cent of the on-pulse data fall within the off- 
pulse gate. However, any monotonic drift in the 
clock, without recovery, would have smeared out 
the pulse obtained by folding the self-data; on 
the contrary we see the on-pulse detected with a 
signal- to-noise ratio of 50-100 (figure[T]). Alterna- 
tively, the clock could have maintained an accurate 
long term average but with large excursions in the 
values of the individual time stamps, i.e. excur- 
sions of the order of 1-2 s (half the pulsar period) 
for an output data rate of one visibility per 0.131s 
or 0.262s. This sounds somewhat contrived; if this 
were the case we should have expected a histogram 
of the time difference between adjacent data to 
show a large scatter. The histograms plotted in 
figure [5] show a tight scatter (rms 0.007s) around 
the expected mean values (0.262s and 0.131s). 

Additionally, any artifact due to time-smearing 
should have resulted in the shorter-period and 
faster-sampled B2045-16 having a higher off-pulse 
fraction, whereas the ratio is much higher for 
B0525+21. 



To determine the probability of finding an unre- 
lated source coinciding with the pulsar we u sed the 
VLA FIRST catalogue (|Becker et al. 19851 ) to de- 
termine the following relationship between source 
counts (N) and flux-density (S): 

logN = 2.2 - 0.826 x logS + log{5S) + log(A sky ). 

(1) 

where S is the flux at 1.4 GHz in mJy, 8S is the 
flux bin width at 1.4 GHz in mJy and A s k y is the 
area of sky under consideration in degree 2 . 

The area of the synthesized beam is ^100 
arcsec 2 at 325 MHz. The flux-density for consid- 
eration in the equation is between 3ximage-rms 
and 5-10 per cent of the pulsar flux, translated to 
1.4 GHz (using spectral index 7 = 0.75). The up- 
per value is related to a realistic upper limit to the 
flux fraction i n the off-pulse from the single-dish 
pulse profile (|Gould fc Lvne 19981 ). This yielded 
a probability of coincident unrelated sources of 
r~> 1.5 x 10~ 3 . The probability of detecting coinci- 
dent sources for both pulsars is the square of the 
above number and therefore highly improbable. 



4.3. Leakage of the signal along the time 
series 

This scenario requires temporal correlation be- 
tween data separated by at least half the pul- 
sar period. The basic sampling interval of the 
GMRT correlator is 31.25 nanoseconds. A cross- 
spectrum is output every 16/is by fourier trans- 
forming 512 such samples. The data output every 
16/j.s should be independent from all other such 
data in the time series. 8188 (or 16376) such inde- 
pendent samples are averaged to obtain the short- 
term acquisition (STA) date output at 0.131 (or 
0.262 s). Given the independence at the 16/j,s level 
it is difficult to envisage correlation between one 
STA and the next, even more so across 6-10 STAs 
(half the pulsar period). Nevertheless this was a 
possibility. Observers using the standard inter- 
ferometric mode of the GMRT would not be af- 
fected as none of their targets vary on sub-second 
timescales, and their integration time is rarely less 
than 2s (more typically 17s). The leakage would 
only redistribute the flux density along the time 
axis without causing any change in the measured 
flux density. Observers using the pulsar mode 
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Fig. 5. — Histograms of the time difference be- 
tween adjacent time-stamps for PSRs B0525+21 
(top) and B2045-16 (bottom). The rms of the 
plots is 7 ms for both. 



would not have detected such an effect because 
they are only sensitive to the relative flux in ex- 
cess of the base level. 

We tested this possibility by estimating the 
temporal correlation in the 0.131s data series. As- 
suming that the off-pulse was entirely due to leak- 
ing of the signal from the on-pulse bin to the off- 
pulse bins we estimated the required leakage for 
a spurious detection of the observed level. There 
was no way of introducing a narrow, clean pulse 
into the GMRT receiver system to directly mea- 
sure such a leakage. So, we investigated the con- 
sequence of such a correlation on the noise data. 

First, we ran simulations to estimate the auto- 
correlation function for noise in the presence of 
leakage sufficient to generate the observed off- 
pulse, i.e. we generated a noise series and smeared 
each data point into subsequent data points ac- 
cording to a particular time profile — e.g. con- 
stant leakage into the next 13 bins; or linearly de- 
creasing leakage into the next 13 bins. The value 
13 is related to the period of the two pulsars which 
covered 14 time bins. 

We then recorded 6x10 minute scans of noise 
data (0.131s integration) but with Front-End ter- 
mination which sealed the telescope at the antenna 
feeds. This ensured that no temporally continuous 
external source (cosmic sources and RFI), which 
would be correlated across the entire observing 
session, was present in the noise data. 

The auto-correlation for each baseline-channel 
(435x120) data was separately calculated for each 
scan. The median autocorrelation profile and the 
scatter are shown in Figure |6l The worst data 
are at the 0.1 per cent level but the median is 
only about 0.04 per cent. The profile does not 
fall to zero at large time lags, indicating non- 
stationarity of the noise signal. We obtained a 
similar behaviour by introducing a time-varying 
mean level into the noise in our simulations. This 
time variation in the mean level may reflect the 
system gain variations as well, which would have 
been corrected to a greater or lesser extent by self- 
calibration during imaging. Therefore the profile 
is a firm upper limit to the contribution of tempo- 
ral leaking of signals to an off-pulse detection. 

It should be noted that interferometric imag- 
ing is an excellent filter of bad data; a compact 
source (like the off-pulse detection) would require 



9 



most of the baselines to show high temporal auto- 
correlation; a few bad baselines with high auto cor- 
relation would only result in noisy ripples across 
the image and not a localised source. Therefore 
the median line is a more accurate measure of the 
leakage for our study than the extreme points in 
the scatter. 

We note that non-stationarity and gain fluctu- 
ation only occur at the 0.05 per cent level in the 
GMRT, and would not have a discernible effect for 
most imaging exercises. 

Figure [7] shows a magnified representation of 
the relevant part of the noise auto-correlation 
plot. Al and A2 (both PSR B0525+21) and Bl 
(PSR B2045-16) indicate the required noise auto- 
correlation if the detected off-pulse were due to the 
leaking of the signal from the on-pulse bin. The 
analysis recipe for PSR B0525+21 (with similar 
arguments for PSR B2045-16) is as follows: 
1. The pulsar period includes 14 time bins, each 
of 0.262 s. The on- and off-pulse fluxes mea- 
sured over 5 time bins are 80.2 and 3.9 mJy, re- 
spectively. Higher t ime resolution observations 



(jGould fc Lvne 19981) show that the main-pulse 
narrower than one of our bins. 

2. A leakage of 0.9 per cent of the pulse flux from 
the first bin into each of the next 13 bins will 
explain the observed off-pulse. This corresponds 
to an auto correlation level of 0.0087. The auto- 
correlation profile calculated for GMRT noise (the 
lowest curve 'N' in figure [7]), is smaller than the 
required amount by a factor of at least 20. 

3. A leakage reducing with temporal distance is 
more likely. A linearly reducing leakage (curve A2) 
requires a similar average auto-correlation in the 
off-pulse region and a higher correlation at short 
lags. 

4. We conclude that temporal leaking of the sig- 
nal is incapable of explaining the detected level of 
off-pulse emission. 

The curve B2 represents the corresponding 
model for PSR B2045-16: off-pulse 4.3 mJy; on- 
pulse flux 305.2 mJy; required constant leakage 
0.43 per cent; required correlation in the off-pulse 
region 0.003; discrepancy is a factor of ~ 6. 

Additionally, the argument against time-smearing 
(Sec. 4.2, second para) is also valid against leaking 
of the signal along the time series. 
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Fig. 6. — Auto correlation profile for front-end 
terminated GMRT noise data. The scatter plot is 
for each baseline-channel-scan while the black line 
is the median profile from the scatter. The non- 
zero value even at large lags is indicative of non- 
stationary noise signal and system gain variation. 
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Fig. 7. — An expanded view of the short-lag re- 
gions of the auto-correlation profile plotted in Fig- 
ure E The plots Al, A2 (PSR B0525+21) and 
Bl (PSR B2045-16) are the required (simulated) 
noise auto-correlation for a spurious off-pulse de- 
tection due to temporal leaking of the pulse into 
subsequent bins. The thick black curve is the 
observed noise auto-correlation profile (from Fig- 
ure [6]). 
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5. Discussion 

5.1. Previous Searches for Off-Pulse Emis- 
sion, including PWNe 

A literature survey yielded many reports of 
searches for "continuous emission in the direction 
of pulsars" (Table |3|) . Here we shall discuss only 
those in which gating techniques were applied to 
target off-pulse emission, which could have been 
either from the pulsar itself or from PWNe. 

Perry k Lyne (1985) reported off-pulse emis- 
sion from 4 pulsars (PSRs B1541+09, B1604- 
00, B1929+10 and B2016+28) using a gated 
two-element interferometer at 408 MHz. Sub- 
sequent studies showed that two of these pul- 
sars, B1541+09 and B1229+10, were aligned ro- 
tators, undermining the claim of off-pulse origin 
( Hankins et al. 1993L Rathnasree k Rankin 1995h . 
The off-pulse emission in the other two pulsars, 
B1604-00 and B2016+28, were later identified 
with unrelated background sources (Storm k Van 
Someren Greve 1990; Hankins et al. 1993). 

Bartel et al. (1985) reported a non-detection 
from PSR B0329+54 using gated VLBI interfer- 
ometry at 2.3 GHz. Strom k Van Someren Greve 
(1990) reported non-detections in 43 pulsars using 
the WSRT in gated mode at 327 and 609 MHz. 
Stappers et al. (1999) reported non-detection in 4 
pulsars using ATCA in gated intcrfcrometric mode 
at 1.3 and 2.2 GHz. 

Gaensler et al. (2000) searched for unpulsed 
emission from 27 energetic and/or high velocity 
pulsars using gating interferometry with the VLA 
and ATCA at 1.4 GHz and found emission in two 
cases. However they could not conclusively de- 
termine whether these detections corresponded to 
PWNe or the pulsar magnetosphere; they even 
suggested that the detections could be spurious 
and a result of improper online gating. These de- 
tections were at the level of 2er (PSR B1634-45) 
and 13ct (PSR B1706-16). 

There have been detections of several ra- 
dio pulsars with associ ated plerionic bow shoc k 
nebula: B 1951+32 dHester k Kulkarni 19881) . 
B1757-24 dFrail k Kulkarni 199ll). B1853+01 



([Frail et al. 19961) and B1643-43 ([Giacani et al. 20011 ) 
All of these are highly energetic pulsars associated 
with SNR (|Chatteriee k Cordes 2o"oi ). Gaensler 
et al. (1998) reported the detection of off-pulse 



emission from PSR B0906-49, which appears to 
be a PWN generated by a slow moving pulsar in 
the dense ISM. 

5.2. Our detections in the context of pre- 
vious effort 

We discuss here the unsuccessful search for off- 
pulse emission in PSR B0525+21 by Weiler et al. 
(1974) and Perry & Lyne (1985). Weiler et al. 
(1974) were looking for extended emission around 
the pulsar without employing any kind of gat- 
ing. They had a synthesized beam size of 47" x 
133"with a noise rms of 2.5 mJy/b. It is appar- 
ent that their setup was inadequate for the present 
level of detection. Perry k Lyne (1985) employed 
a two-element interferometer with hardware gat- 
ing at 408 MHz for their studies but could neither 
image the fields and nor correct for sensitivity vari- 
ations. They listed considerable emission in the 
off-pulse bin for B0525+21 (in fact greater than 
for B1929+10 and B2016+28, which were claimed 
to be positive detections) but the noise was also 
very high precluding a positive detection. 

The major advantages of our effort were: 

1. Low frequency of observation: the steep spec- 
trum pulsar is likely to dominate over emission 
from any associated nebula with flatter spectrum; 

2. A multi-element interferometer: The GMRT 
is currently the most sensitive low-frequency in- 
strument. The ability to correct gain variations, 
filter bad dat a using imagin g residuals and RFI 
excision tools (lAthreva 2009h allowed us to make 
sufficiently deep images. 

3. Off-line gating of data: this allowed secure sep- 
aration of the off- and on-pulse sections. 

5.3. PWN emission 

The confinement of the relativistic wind from 
pulsars generates PWNe which are luminous 
across the electromagnetic spectrum in syn- 
chrotron, inverse compton, and optical line emis- 
sion from the shocked regions. All known PWNe 
seen around radio pulsars have spin-down lumi- 
nosities E > 10 34 erg s _1 . They all appear to be 
young and, with the exception of PSR B0906-49 
(jGaensler et al. 19981 ). are associated with SNR. 
Two classes of PWNe can be formed in the ab- 
sence of associated SNR: one, where the pulsar 
wind is confined by the density of the ISM - 
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Table 3: Literature survey for off-pulse studies. 



Non detections/ 
subsequent SNR 


Detections 
later refuted 


Tentative 
detection 


PWNc detection 


1. Schonhardt 1973, 1974 

2. Wciler ot al. 1974 

3. Cohen ct al. 1983 

4. Bartel et al. 1985 (G) 


1. Gopal-Krishna 1978 

2. Glushak ot al. 1981 

3. Perry & Lyne 1985 (G) 


1. Gacnslcr ct al. 2000 (G) 


1. Frail & Kulkarni 1991 

2. Frail et al. 1996 

3. Gacnslcr ct al. 1998 

4. Giacani ct al. 2001 



5. Strom & Van Somcrcn 
Greve 1990 (G) 

6. Hankins ct al. 1993 (G) 

7. Frail & 

Scharringhauscn 1997 

8. Stappcrs ct al. 1999 (G) 



There have been no unambiguous detection of off-pulse emission not associated with PWNe. The gated interfcrometry studies 
in table [3] have been indicated with G. 



known as static PWN; and the other, where the 
wind is confined by the ram-pressure of motion of 
the pulsar through the ISM — called bow-shock 
nebula. Our two target pulsars, B0525+21 and 
B2045-16, are not associated with SNRs. So, we 
explored the possibility that the detected off-pulse 
emission was due to a PWN. 

The efficiency factor ep is defined as the ratio of 
the radio bolometric luminosity (Ln) of a PWN to 
the spin-down energy (E) of the pulsar Lr = epE. 
If we assume a typical PWN spectral index of 7 ~ 
0.3 its radio luminosity between 10 7 Hz and 10 11 
Hz is given by L R = 3.06 x 10 28 d 2 kpc S mJy erg s -1 , 
where dkpc is the distance to the PWN in kpc 
and S m jy is the integrated flux of PWN at 325 
MHz in mJy. Using the measured flux-density 
(table [5]) , source distance and the spin-down en- 
ergies (table [lj we calculate efficiency factors of 
2xl0~ 2 for B0525+21 and 2xl0~ 3 for B2045-16. 
These €r values are 1-2 orders of magnitude higher 
than those of previously known PWNe (typical 
6^ ~ 10~ 4 ; see Frail & Scharringhausen 1997, 
Gaensler et al. 2000 for a discussion). 

Static PWN The relativistic particles and 
Poynting flux emanating from the pulsar, at rest 
relative to the ISM, will drive through the am- 
bient mediu m a shock of radius R s given by 
(|Arons 19831) 



The velocity of the shock front is given by 

R S = U J,) V5 = 3.3 ( lh km s- 

(3) 

Here p = m#n, where m# is the proton mass and 
n is the particle density of the ISM. Using equation 
@ the required particle density for a PWN is 

n = 5.35x 10 11 [ j cm" 3 (4) 

\ ^0.01 / 

where E31 is the spindown power in units of 
10 31 erg s _1 , te is the age in units of 10 6 yr, 
i?o.oi is the radius of PWN in units of 0.01 pc and 
no. 01 is the ISM density in units of 0.01 cm -3 . In 
the present exercise we assume that the PWNe in 
our targets are smaller than one synthesized beam 
width since the off-pulse emission is unresolved in 
both pulsars. Using equation ((4]) and values in 
tables Q] and [5] we determined the ISM density re- 
quired to drive a PWN to be ~ 3 x 10 9 cm -3 for 
B0525+21 and - 10 12 cm" 3 for B2045-16. The 
required particle densities are much higher than 
typical ISM densities of ~ 0.03 cm~ 3 , suggesting 
that these pulsars are too weak and old to power a 
static PWN through the ISM. Since our estimate 
of the sizes of the nebulae are upper limits the 
particle density estimates are lower limits. 

Bow Shock PWN This requires that the shock 
velocity (R s ) be much smaller than the pulsar 
transverse velocity (Vpsr)- The typical ISM den- 
sity of 0.03cm -3 in equation yields R s ~ 



12 



2.8 km s- 1 for B0525+21 and ~ 2.5 km s" 1 for 
B2045-16, which are indeed much smaller than the 
Vpsr list ed in Table ffl The radius of th e shock is 
given by (jFrail fe Scharringhausen 1997D 



RBi 



E 



1/2 



^CPoVpSR 



(5) 



From equation (|5|) the particle density required 
to sustain a bow shock PWN is 



n = 1.67 x 10" 



E- 



31 



R 2 V 2 
n-o.Ol v loo 



(0) 



where E31 is the spin-down power in units of 
10 31 erg s _1 , i?o.oi is the radius of PWN in units 
of 0.01 pc and V100 is the velocity of the pulsar 
through the ISM in units of 100 km s -1 . This 
requires an ISM density of ~ 5 x 10~ 6 cm -3 for 
B0525+21 and - 8 x 10~ 6 enr 3 for B2045-16. 
These values are 3-4 orders of magnitude lower 
than the typical ISM density. This implies that the 
pulsars are too weak to drive a bow shock nebula 
of the size corresponding to the telescope resolu- 
tion. However, a bow shock interpretation can be 
salvaged if the nebulae are ~2 orders of magnitude 
smaller. 

In summary, identifying the off-pulse emission 
detected here with PWNe results in unrealistic val- 
ues of the ISM particle density. This is not sur- 
prising because we selected as our targets middle- 
aged pulsars, which category has hitherto not been 
known to be associated with PWNe. 

5.4. Emission from the Magnetosphere 

Numerous published studies provide estimates 
of a and for the pulsars studied here. Some use 
an empirical/geometrical (E/G) approach to es- 
tablish a and (Lyne & Manchester 1988; Rankin 
1993; see Everett & Weisberg 2001 for a discussion 
on the E/G approach). One can also obtain a and 
values by fitting the RVM to the PPA traverse 
(jMitra fe Li 2004h . but these fits yield highly cor- 
related estimates of the two angles. This issue has 
been discussed in depth by Everett & Weisberg 
(2001) and Mitra & Li (2004). 

Everett & Weisberg (2001) list six different 
studies of PSR B0525+21 with a ranging from 



116° to 159°lJ Mitra & Li (2004) estimated 
a ~ 127 - 144° for PSR B2045-16. Both RVM 
fit and the E/G approach were used to obtain the 
above estimates of a. These studies provide good 
evidence that our targets, B0525+21 and B2045- 
16, are not aligned rotators and therefore any off- 
pulse emission must arise far from the magnetic 
pole. 

The pulse profiles of both these pulsars show 
sharp rising and falling edges. PSR B0525+21 
is an example of a classical double (D) profile, 
whereas PSR B2045-16 is classified as a triple 
([Rankin 1 993). The main-pulse width for these 
pulsars span about 5 per cent of the period, and 
there are no obvious emission components visible 
outside the main-pulse. Based on the pulse widths 
and the estimates of a and 0, Mitra & Li (2004) 
concluded that the radio emission arises at a dis- 
tance of 1-2 per cent of the light cylinder. Hence, 
the radio on-pulse emission of PSR B0525+21 and 
B204516 are classical examples of emission from 
open dipolar field lines. The remaining 95 per 
cent of the period, i.e. the off-pulse region, were 
hitherto thought to be radio quiet zones of the 
pulsar magnetosphere dGoldreich fc Julian 19691 



Ruderman fc Sutherland 19751 ). 

Our off-pulse phase corresponds to 80-208° (from 
the peak of the main-pulse) for B0525+21 and 101- 
229°for B2045-16. In a few cases pre-/post-cursor 
emissions have been detected about 60°from the 
main-pulse, posing ques tions regarding their origin 
(jMitra fc Rankin 2 004). Our detections are even 
further away from the main-pulse. It remains to 
be seen if these off-pulse detections are an ex- 
treme example of PPC or if it represents impulsed 
emission from throughout the pulse period. 

It is to be noted that the opening angle of dipo- 
lar field lines scales as the square of the emission 
height. If radio emission from pulsars originate 
from higher up the light cylinder (3> 1-2 per cent of 
the light cylinder), the polar emission could span a 
larger longitudinal range. This implies that radio 
emission at larger distances from the neutron star 
surface are potential sites of PPC or off-pulse emis- 
sion in pulsars. Alternatively, if the emission arises 



following Everett &; Weisberg (2001) a values should span 
0-180° (and not 0-90°) for a consistent definition of the pul- 
sar rotation axis. Everett & Weisberg provide a compila- 
tion of published a and /3 values corrected for this conven- 
tion. 
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close to the neutron star surface, the emitting re- 
gion will need to encroach upon the closed field 
lines of the neutron star. Thus, the existence of 
magnetospheric off-pulse emission should prove to 
be an important input for pulsar electrodynamic 
models which try to establish the relation between 
magnetos pheric currents an d coherent pulsar radio 
emission (jsnitkovskv 2006h . 



6. Conclusion 

We report GMRT observations at 325 MHz 
which detected off-pulse emission from two long- 
period pulsars which have low spin-down rates and 
are not associated with supernova remnants. We 
have adduced evidence that the signals are not ar- 
tifacts of the observing procedure, nor unrelated 
background sources which happen to be coinci- 
dent with the pulsar. We have argued that ex- 
plaining these detections as Pulsar Wind Nebula 
requires ISM particle densities which differ from 
typical measured values by several orders of mag- 
nitude. Robust estimates of the geometrical pa- 
rameters of these two pulsars argue against them 
being aligned rotators. This leaves the possibility 
of emission from the magnetosphere. If the off- 
pulse emission arises from a much higher height 
than the on-pulse then it could arise from open 
field lines at the edge of the light cylinder. On 
the other hand if the off-pulse emission is from 
the same height as the on-pulse emission then it 
must be associated with closed field lines. Fur- 
ther studies at multiple frequencies, in polariza- 
tion and higher resolution are needed to establish 
the nature of this emission which can impose valu- 
able constraints on the pulsar emission mechanism 
models. 

We thank the director of NCRA for discre- 
tionary observing time which made this work pos- 
sible. We also thank the staff of the GMRT for 
their help with these non-standard observations 
possible. The GMRT is run by the National Cen- 
tre for Radio Astrophysics of the Tata Institute of 
Fundamental Research. We also thank the referee, 
Scott Ransom, for useful comments. 
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